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ABSTRACT 

We report the discovery of two simultaneous quasi-periodic oscillations in the per- 
sistent flux of GX 349+2 at frequencies 712 ± 9 and 978 ± 9 Hz, with rms amplitudes 
1.25%±0.34% and 1.34±0.32%, respectively. During our 152 ks observation with the 
Rossi X-ray Timing Explorer, GX 342+2 was in either the normal branch or the flaring 
branch with count rates in the nominal 2-60 keV RXTE-PCA band ranging from a low 
of 8,000 cps to a high of 15,000 cps. The kHz QPOs were observed only when the source 
was at the top of the normal branch when the count rate was about 8,200 cps corre- 
sponding to a flux of 1.4 x 10~ 8 ergs cm -1 s -1 in the 2-10 keV band. With this report, 
now kHz QPOs have been observed in all the 6 Z sources. 

Subject headings: accretion, accretion disks — stars: neutron — stars:individual (GX 
349+2)— X-rays:stars 
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1. Introduction 

GX 349+2, also known as Sco X-2, is one of 6 
bright low mass X-ray binaries whose X-ray color- 
color diagrams resemble the shape of the letter Z, 
which, t herefore, have been referred t o as the Z 
sources ( Hasingcr fc van der Klis 1989 ). They are 



among the brightest X-ray sources in the sky. Their 
time variability in the frequency domain below 100 Hz 
as charatcerized by their FFT power spectra is closely 
correlated with their energy spectral state as charac- 
terized by their X-ray color-color diagram. Since the 
launch of the Rossi X-Ray Timing Explorer (RXTE) , 
quasi-periodic intensity oscillations with frequencies 
in the range of several hundred hertz to over one 
thousand hertz have been discovered in 5 of the 6 



Z sources: Sco X-l (van der Klis et al. 1996), Cyg 



X-2 (Wijnands et al. 1998), GX 17+2 (Wijnands et 



al. 19971), GX 5-1 Jva,n der Klis et al. 199G| ), and GX 
340+0 jwijnands et al. 199S| ). GX 349+2 is the only 
source from which no kHz QPOs have been reported 
to date. 

In this letter we report the discovery of kHz QPOs 
in GX 349+2 and compare their characteristics with 
those of the QPOs observed in the other Z and atoll 
sources. 

2. Observations 

The Rossi X-Ray Timing Explorer (RXTE) was 
launched on 30 December 1995. It carries three X- 
ray instruments: an array of 5 xenon gas proportional 
counters that has a nominal bandwidth of 2-60 keV, 
8 Nal and Csl scintillation detectors sensitive to X- 
rays in the 15-250 keV band, and an all sky monitor 
(ASM) sensitive to X-rays in the 1.5-10 keV band. 
The sailent features of the RXTE satellite are its 
large X-ray collection area (6,250 cm 2 for the PCA 
alone) and its large telemetry beandwidth enabling 
time-tagging each X-ray to a micro-second accuracy. 
For the results reported in this paper, we use only the 
proportional counter array (PCA). 

The observation took place between 9 and 29 Jan- 
uary 1998 with a total accumulated exposure time 
of 152 ks. It was broken up into 55 pieces by the 
satellite observation planning process, the South At- 
lantic Anomaly, and Earth occult of the source, with 
the shortest being only 400 s and longest 3,900s, with 
2,800s being the typical duration. 

In addition to the two standard data modes that 



are available to all RXTE-PCA observations, we used 
four additional data modes: three single-bit modes 
each with a time resolution of 128/it (lfxt = 2~ 20 s) 
and cover the PCA pulse height channels 1-13, 14-17, 
and 18-23, which correspond to energy bands, respec- 
tively, 2-5.0, 5.0-6.5, and 6.5-9.0 keV. Photons with 
energies above channel 23 (9.0keV) are recorded with 
an event mode which has a time resolution of 128 /it 
and 64 energy channels. 

3. Data Analysis and Results 

We constructed the color-color diagram as shown 
in Figure [l] using the Standard-II data and the stan- 
dard PCA background model. Then we combined the 
3 single-bit data streams and the event mode data 
stream to form a single time series with a time reso- 
lution of 256)it to construct an FFT power spectrum 
for every 32s of the data. We then obtained the aver- 
age power spectrum for each of the 55 segments and 
rebinned it from a frequency resolution of 1 /32 Hz by 
a factor of 512 to 16 Hz per bin. Visually scanning all 
the 55 average powered spectra, we found one of the 
55 segments showed two significant peaks at 706 Hz 
and 999 Hz, respectively. All the other segments did 
not have any significant peaks anywhere above 100 
Hz. 

We then tried the following procedures to enhance 
the detection significance of the peaks: (1) summing 
up all the 55 segments, (2) summing up only those 
segments with count rates below 9,000 cps, (3) re- 
moving events in the first single-bit data stream, i.e., 
those events with pulse height channel below 14, and 
(4) summing up those segments at the top of the nor- 
mal branch in the color-color diagram. Procedures 
(1) and (2) washed out the peaks completely. Pro- 
cedure (3) did not measureably alter the final signifi- 
cance. Only procedure (4) significantly enhanced the 
significance. The resulting power spectrum is shown 
in Figure g. 

We fitted the power spectrum in Figure || to two 
Gaussian peaks plus a Poisson noise term: 
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The best fit parameters are listed in Table [[]. Their 
errors correspond to a change A% 2 = 1. 

The rms amplitude is calculated according to the 
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following formula: 



i^Ao 65536 



V N ph 2048 ' (2) 

where A and tr are as defined in equation |l|; N p ^ = 
263, 840 is the number of photons in 32s averaged 
over the entire time interval from which Figure || 
is constructed, and the factor 6,5536/2,048 converts 
the the power amplitudes from per frequency inter- 
val to the original FFT amplitudes ( |van der KUe 
198^|). Plugging in all the numbers, we obtain the the 
rms amplitudes for the two QPOs, 1.21 ± 0.33% and 
1.30 ± 0.31%, where we have folded both the errors 
in A and the errors in a into the final rms amplitude 
errors. 

There are two factors that have systematically 
suppressed the measured rms with respect to the 
true values: detector background events and detec- 
tor deadtimc. Background event rate only consti- 



tutpd no more than 2% of the overall count rate, 



which translates into 1% error on the rms amplitude. 
Given that the PCA detectors has a deadtime of 10/is 



(Zhang 1995) and each detector had a count rate of 



8,000/5 cps, the overall correction to the rms am- 
plitudes amounts approximately to 2% (van der Klis 
198^|). In other words, combining the effects of de- 
tector background events and the detector deadtime, 
we should revise up the above quoted rms amplitudes 
by 3%, i.e., multiply them by a factor 1.03. The final 
rms amplitudes for the two QPOs are: 1.25±0.34% 
and 1.34±0.32%. 

We note that there is an intriguing, but statisti- 
cally insignificant peak at 1,020 Hz in the GX 349+2 
power spectrum reported by Kuulkers and van der 
Klis (1998). With the hindsight gained with our re- 
sults, we think that their 2.6a peak most likely is due 
to the kilohertz QPOs in the flux, not due to a sta- 
tistical fluctuation. 

4. Discussion 

We have reported the discovery of two simultane- 
ous QPO peaks in the FFT power spectrum of GX 
349+2. In this section we place our results in the con- 
text of the kHz QPO phenomenologies of the other Z 
sources. 

GX 349+2 is the last Z source from which such 
characteristic twin QPO peaks have been found. The 
color-color diagram we observed is also quite charac- 
teristic of those observed in the past. Compared to 



other Z sources, GX 349+2 has never been observed 
to be on its horizontal branch. Four of the 6 Z sources, 
i.e., Cyg X-2, GX 5-1, GX 17+2, and GX 340+0, have 
clearly observable horizontal branches. The kHz QPO 
characteristics of those sources are nearly identical in 
that they all show up on the horizontal branch of the 
Z diagram. As the source moves along the Z track to- 
ward the vertex of the horizontal branch and the nor- 
mal branch and further down the normal branch, the 
QPO centroid frequencies increase and their rms am- 
plitudes decreases until they cease to be observable in 
the middle part of the normal branch. The frequency 
difference of the twin QPO peaks are consistent with 
being constant for all of these sources. The differ- 
ences are, respectively, 343+21, 327+11, 294+8, and 
374+24 Hz. In the case of Sco X-l, the kHz QPO 
peaks are observed on the normal branch. The sepa- 
ration of the two QPO peaks changes from a high of 
292+2 to a low of 247+3 Hz as Sco X-l moves down 



the normal branch toward the flaring branch (van dei 



Klis et al. 1997b) 



The horizontal branch, if any, of the Z track of the 
GX 349+2 has not been observed. The kilohertz QPO 
characteristic we have observed in GX 349+2 is very 
similar to those of Sco X-l and the other Z sources 
in that the kHz QPOs exist on the upper or middle 
normal branch. As the source approaches the flaring 
branch, the kilohertz QPOs become weak and cease to 
be observable. Unfortunately, our existing data from 
GX 349+2 are not sufficient in quantity to further 
investigate the characteristics of the observed QPOs. 
Further and more extensive observations, especially 
when GX 349+2 is at the upper normal branch, are 
highly desirable and encouraged. 

We thank E. H. Morgan for helping us setting up 
the PCA data modes. We made use of the software 
and database made available by the HEASARC at the 
Goddard Space Flight Center and the CERN software 
library. 
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Table 1: Parameters of the best fit in Figure 0. The 
errors quoted all correspond to a change of \ 2 by 1. 



Parameters 


Best Estimates and Errors 
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At 


0.0154 ±0.0036 


h (Hz) 


712.2 ±8.8 


o-i (Hz) 
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0.0142 ±0.0030 


h (Hz) 


978.0 ±9.5 


0-2 (Hz) 


39.1 ±9.9 


C 


1.950 ±0.001 
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Fig. 1. — Color-color diagram of the 55 segments of 
data. The filled and open circles represent those seg- 
ments with and without observable kHz QPO peaks 
in their FFT power spectra, respectively. 
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Fig. 2. — FFT power spectrum from the three seg- 
ments of data as indicated by filled circles in Figure [j]. 
The smooth line connecting the points is the best fit 
with two Gaussians plus a constant. See text for the 
best fit parameters. 
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